Diffusion on a heptagonal lattice 
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We study the difFusion phenomena on the negatively curved surface made up of congruent hep- 
tagons. Unhke the usual two-dimensional plane, this structure makes the boundary increase ex- 
ponentially with the distance from the center, and hence the displacement of a classical random 
walker increases linearly in time. The diffusion of a quantum particle put on the heptagonal lattice 
is also studied in the framework of the tight-binding model Hamiltonian, and we again find the 
linear diffusion like the classical random walk. A comparison with diffusion on complex networks is 
also made. 
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The transport phenomena have been extensively stud- 
ied on various spatial structures to reveal the character- 
istic thermal and electrical properties of materials. The 
random walk is one of the classical topics for study- 
ing those phenomena in the nonequilibrium statistical 
physics [l|, and the transport on highly disordered sys- 
tems such as complex networks is also of recent research 
interest 0, B S|- Since the average path length is of- 
ten greatly reduced in complex networks, one can easily 
expect that a particle diffuses there faster than on the 
ordinary regular lattice structures. 

Recently, the surface with a constant negative Gaus- 
sian curvature is under active investigation to study geo- 
metrical effects on the critical phenomena 0, Q ■ A neg- 
ative Gaussian curvature means that the surface locally 
looks like a uniform saddle everywhere, with the resulting 
surface getting highly curled as one moves outward [7|. 
Consequently, the boundary length of a circle with ra- 
dius r is given by 27r sinh r, rather than by 27rr [8|, 
which means in turn that the characteristic path length 
increases only logarithmically with the system size N Q . 
For a positively curved surface, the system becomes even- 
tually closed, like a sphere, where the magnitude of cur- 
vature has to approach zero in the thermodynamic limit. 
Accordingly, there is little reason to study this case sepa- 
rately. In contrast, surfaces of a constant negative Gaus- 
sian curvature can be extended indefinitely, which makes 
such a geometry apt for studying novel physical proper- 
ties. Furthermore, the development of nanotechnology 
can make it possible to construct such a structure in re- 
ality, and the physical properties discovered here are also 
of practical importance. There already exist theoretical 
interests in negatively curved nanostructures with a very 
large surface area embedded in a limited volume [9| . 

It has been reported that the hyperbolic Brownian mo- 
tion has a limited direction [l3| with a constant outward 
drift [ll| . Such nontrivial behaviors are nothing myste- 
rious and directly related to the exponentially increasing 
length of a boundary: Since there is always a larger num- 
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FIG. 1: (a) Heptagonal lattice represented on the Poincare 
disk with the total number of levels L — 4. The metric is 
given in such a way that all heptagons are congruent and the 
circle indicates the points at infinity, (b) Random walk on the 
heptagonal lattice where a walker at the black dot has three 
choices (colored in gray) at each time, (c) In the dual lattice, 
walking around a heptagon is equivalent to the reflection of a 
triangle. 



ber of points outside than inside by some constant pro- 
portion, a random walker tends to move outward, which 
allows the random walker only a limited direction 12]. 
However, the validity of such a simple geometric under- 
standing of the fast classical diffusion needs to be taken 
carefully if one compares the classical diffusion with the 
quantum one, which composes the main motivation of 
this Brief Report. 

In this work, we present how to simulate the random 
walk on such a geometry and reconfirm the linear classical 
diffusion for a heptagonal lattice structure. The quantum 
diffusion of a particle is then investigated in compari- 
son both with the classical diffusion behavior and with 
Ref. [4] in which the quantum difFusion time has been 
shown to scale as Tq ~ N"^^ while Tc ^ N" for the 
classical diffusion. Also in Ref. [3|, quantum mechanical 
transport on graphs has been shown to be faster than 
the classical one, except on some finite treelike graphs. 
In comparison, Tq ~ Tc is revealed in the present work on 
the heptagonal lattice geometry. It is remarkable that 
one can drastically accelerate the transport by introduc- 
ing heptagonal plaquettes to the lattice [9|. 

Let us describe the structure of the heptagonal lattice. 
Similarly to the two-dimensional (2D) plane which can be 
covered by congruent triangles, squares, and hexagons, 
a negatively curved surface can be tiled by polygons 



through the hyperbohc tessellation [l3|: Suppose a lat- 
tice where q regular p-gons meet at each vertex. For 
example, the triangular lattice is denoted as p = 3 and 
g = 6. It is known that the tiling with given p and q 
covers the negatively curved surface if {p — 2) (g — 2) > 4. 
If we choose its basic element as the simplest polygon, 
the regular triangle, at least seven triangles should meet 
at each vertex leading to {p, q\ = {3, 7}. Taking its dual 
lattice, we obtain the heptagonal lattice structure with 
{Pi ?} = {7, 3}. The resulting structure can be most suit- 
ably represented in the Poincare disk [14| as in Fig.[IJa). 
Since the heptagons are located in a concentric fashion, 
we denote the central heptagon as the first level. Accord- 
ingly, its seven nearest-neighbor heptagons constitute the 
second level, and the next nearest ones the third level, 
and so on. Figure [IJa) thus describes the heptagonal 
lattice up to the level I = 4. Let H{1) be the set of hep- 
tagons in the Zth level. The set of vertices on the outward 
boundary of H{1) is called the lih layer. As one can see, 
this structure uniformly and completely fills the surface 
and provides an adequate tool for studying physics of 
the negative Gaussian curvature. For numerical calcula- 
tions, we may use the heptagonal lattice constructed up 
to some finite level, L, and capture the finite size effects 
by varying L. 

We first study numerically the random walk problem 
in our heptagonal lattice structure, and the results are 
compared with numerical solutions of the diffusion equa- 
tion for the hyperbolic geometry. A possible way of sim- 
ulating a random walker on the heptagonal lattice is to 
identify the connection structure first and then let a par- 
ticle move along the edges. In Fig. EJJa), the position 
of the random walker measured by the level I is shown 
as a function of the time step n for the heptagonal lat- 
tices of sizes L = 7, 10, and 14. It is clearly shown that 
the random walker drifts away from the starting position 
I = 1 linearly in time. As the lattice becomes larger (as 
L is increased), the linear diffusion regime becomes ex- 
tended, indicating (/) ex n in the thermodynamic limit. 
One drawback of this approach is that the number of 
points increases exponentially as we add concentric lay- 
ers. The memory constraint restricts the distance from 
the origin, and therefore the results severely suffer from 
finiteness of the model system under consideration. 

A better alternative is obtained from the hyperbolic 
tessellation. Since the heptagonal lattice is dual to {3, 7}, 
the transition between the neighboring points coincides 
with reflecting a hyperbohc triangle [see Fig.[TJc)]. Since 
the Poincare disk can be identified with a unit disk on the 
complex plane, we start from a regular triangle (zi , Z2, 2:3) 
where Zi's are complex numbers. Let one of the vertices, 
say zi , be located at the origin. Setting the interior angle 
at that point to be 27r/7, the triangle is regular only when 
\zi — Z2I = \zi — Z3I « 0.49697, as this surface has its own 
intrinsic length scale, i.e., the curvature [IJI- If a triangle 
[zi, Zj, Zk) is reflected around the edge along Zi and Zj, we 
get a new triangle (z^, Zj, z'f.) with z'f. = w ■ 
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FIG. 2: (Color online) (a) Average level (/) of the random 
walker versus the time n. The heptagonal lattice of the size 
L is first built, and we perform random walks as in Fig. [Ub). 
We see (Z) oc n for small n, while (l) eventually saturates 
as the walker approaches the external boundary, (b) Aver- 
age distance {dh) from the origin versus the time n, obtained 
from the method in Fig. [TJc). The deviation from the linear 
behavior is not due to the finiteness of the lattice as in (a), 
but due to accumulated numerical errors. All the averages 
were taken from 10^ independent realizations. 



^ = jw — Zi| = |i(j — Zjl, and z is the complex conju- 
gate of z. The center of a triangle, z, has the hyperbolic 
distance from the origin by dh{z) = ln[(H- |z|)/(l— \z\)]. 
Note that dh diverges to infinity as \z\ —> 1, which is 
consistent with the definition of the Poincare disk. The 
result is depicted in Fig. [2Kb), which clearly shows that 
the distance is linearly proportional to time. The numer- 
ical inaccuracy becomes larger as we repeat the reflec- 
tion of triangles. In order to confirm the origin of the 
deviation from the linear behavior (dh) ^ n at large n, 
we intentionally assign numerical precision and observe 
how {dh) depends on it. In Fig. [2Jb), it is seen that as 
we use better numerical precision, the linear regime be- 
comes more extended, indicating that the deviation from 
the linear behavior originates from the simple artifact of 
the numerical accuracy. Beyond the deviation point of 
each curve in Fig. ^ the numerical values are not trust- 
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FIG. 3: (Color online) The solution (/!)(r, t) of the hyperbolic 
diffusion equation. Inset: Average distance from the origin 
(r) = J rcj>{r-) sinh r dr increases linearly in time. 



worthy and thus the decrease of {dh) should not be taken 
as real. 

For the random walk on a 2D flat surface, the walker 
has zero probabiUty of escape by Polya's theorem which 
states that the walk becomes transient for dimensions 
larger than two [l] . If we denote the probability density of 
flnding a particle at a position r after time t as 0(r, i), the 
diffusion process is described by d(j)/dt — V^0. Solving 
this equation by Fourier transformation, one can see that 
the expected displacement from the origin is proportional 
to the square root of time in the 2D flat surface. The 
Laplace operator V^ is changed to the Laplace-Beltrami 
operator for the hyperbolic metric [15J . 
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and we get the solution of the hyperbolic diffusion equa- 
tion, d(j)/dt^ A(j) [ll[3: 
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with the distance p between two positions (r, &) and 

(r', B'), and /(i, p) = /""(se-^'/^Wcosh s 



cosh p)ds. 

The numerical integration is then performed to get the 
result presented in Fig. [31 which again shows the linear 
diffusion (r) ex t. We thus conclude that our different 
approaches, i.e., discrete random walks by triangle re- 
flections on the Poincare disk and the solution of the 
hyperbolic diffusion equation, unanimously confirm that 
the distance from the origin increases linearly in time. In 
comparison to (d) ~ \ft for the 2D flat surface, the diffu- 
sion occurs much faster in the negatively curved surface. 
We next study the diffusion of a tight-binding quantum 
particle in the heptagonal lattice. The quantum diffusion 
phenomena, and the existence and nature of the localiza- 
tion transition have been studied on regular lattice struc- 
tures for a long time. Motivated by the intensive research 



FIG. 4: (Color online) Diffusion of the tight-binding quantum 
particle on the heptagonal lattice. As the time t evolves, the 
quantum particle diffuses toward the upper layers. Note that 
the linear regime extends as the size becomes larger. 



interest of complex networks, recent years have observed 
the beginning of research on quantum mechanical sys - 
tems put on the structure of complex networks [3, l3|, |l7| . 
We assume that the wave function is localized on the 
first layer containing seven points [see Fig. [Ha)]. The 
time evolution of the quantum particle is governed by 
the Schrodinger equation i{d/dt)\ilj) — 'H\4') {h = 1). By 
using the perfectly localized states as basis kets, we ap- 
ply the tight-binding approximation that Tiij = 1(0) if 
two points i and j are connected (not connected). The 
Schrodinger equation is then numerically integrated by 
the fourth order Runge-Kutta method. The spreading 
of the wave packet is measured by the average layer, (/), 
depicted in Fig. [4] Again, one can see clearly that the 
particle diffuses linearly in time as in the classical dif- 
fusion. The saturation comes from the finiteness of our 
lattice. It is to be noted that in the negatively curved 
heptagonal lattice, both the quantum and the classical 
diffusion exhibit the linear diffusion. 

We then take an alternative approach to investigate 
the quantum diffusion problem, through the use of the 
mapping of the above tight-binding Hamiltonian to the 
one for a free particle. We note that the shift of energy 
by a constant amount does not change any measurable 
quantity and makes the transformation Ti. -^ Ti — kl with 
a degree k (the number of neighbors, e.g., k — 3 for the 
heptagonal lattice) and the identity operator /. This sim- 
ple transformation makes the Hamiltonian proportional 
to the lattice Laplacian, except for the outermost points 
where fc 7^ 3. In short, the tight-binding Hamiltonian 
can be phenomenologically treated as that of the free 
particle simulated on the discrete lattice, described by 
id'^/dt = Avl" with the wave function \E' = (r|^). Ac- 
cordingly, if we substitute the time t in the Schrodinger 
equation by the imaginary time it, it takes exactly the 
same form as the diffusion equation. It is then straight- 
forward to apply the previous result on classical diffusion 




FIG. 5: (Color online) Solution of the quantum free particle 
Schrodinger equation in the negatively curved surface. Inset: 
Average distance from the origin (r) = J r\'^{r)\'^ sinh r dr is 
shown to increase linearly in time. 



in hyperbolic lattice (see Fig. [3]) to obtain the propagat- 
ing solution of the quantum particle on the hyperbolic 
plane ^16i] (see Fig. [5]). It is revealed that the average 
distance again exhibits the linear diffusion property as 
shown in the inset of Fig. [5] in accordance with Fig. H) 

In the Euclidean or disordered structures, the quantum 
diffusion has been known to be faster than the classical 
one 0, Q. In contrast, the classical motion now becomes 
comparable to the quantum one in speed on the neg- 
atively curved surface, with the aid of the geometrical 
drift. Some differences, however, seem to remain: The 
geometrical drift pushes the particle only in the outward 
direction from the starting point. Based on the report for 
tree structures, a quantum particle is believed to readily 
propagate back to the origin [l^. In fact, a classically 



diffusing particle acts somewhat like a flow, even if the 
initial direction may be chosen randomly, in the sense 
that it appears to have some nonzero velocity. If one in- 
troduces a real flow here, however, the flux would soon 
become negligible because of the exponentially increasing 
boundary. Our observation is qualitatively similar to the 
report for the Bethe lattice [ll|, since the angular move- 
ment is effectively suppressed in the long run. Yet it is 
still notable that the heptagonal lattice provides a better 
representation for the negatively curved surface, and that 
the numerical simulation could be taken with more ease 
by the hyperbolic tessellation technique explained above. 
One may expect that the transport can be enhanced by 
introducing a negative Gaussian curvature in biological 
or engineering applications. 

In summary, we have investigated the classical and 
quantum diffusions in the heptagonal lattice, which is 
a discrete representation of the surface with a constant 
negative curvature. Even a classical particle has been 
shown to diffuse so fast that the average displacement 
is linearly proportional to time. The quantum diffusion 
has also been shown to exhibit the same linear diffusion 
behavior. Those results on discrete lattices were also con- 
firmed in a continuum by solving a hyperbolic diffusion 
equation in real and imaginary times for classical and 
quantum diffusion, respectively. 
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